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ABSTRACT 

The outer regions of galactic disks have received increased attention since ultraviolet observations 
with GALEX demonstrated that nearly 30% of galaxies have UV emission beyond their optical extents, 
indicating star formation activity. These galaxies have been termed extended UV (XUV) disks. Here, 
we address whether these observations contradict the gas surface density threshold for star formation 
inferred from Ha radial profiles of galaxies. We run smoothed particle hydrodynamics simulations of 
isolated disk galaxies with fiducial star formation prescriptions and show that over-densities owing to 
the presence of spiral structure can induce star formation in extended gas disks. For direct comparison 
with observations, we use the 3-D radiative transfer code Sunrise to create simulated FUV and K s - 
band images. We find that galaxies classified as Type I XUV disks are a natural consequence of spiral 
patterns, but we are unable to reproduce Type II XUV disks. We also compare our results to studies 
of the Kennicutt-Schmidt relation in outer disks. 

Subject headings: galaxies: spiral, galaxies: structure, galaxies: evolution, ultraviolet: galaxies 



1. INTRODUCTION 

Star formation rate profiles derived from Ha emission 
indicate that galaxies follow a "Kennicutt-Schmidt Law" 
where the star formation rate and gas surface densities 
averaged in azimuthal annuli are related by Ssfr ^g& s 
([Kennicuttl Il998bl [l989). However, sudden drops in 
Ha intensity create strong departures at a gas surface 
density of 3-5 M pc -2 . This is generally interpreted 
as a threshold density for star formation ([Kennicuttl 
1989; iMartin fc Kennicuttl l2QQl| ). Several explanations 
for this threshold have been proposed, including a tran- 
sition between dy namically un stable and stable regions 
of the galax y (e.g. iToo mre 196 j) or a phase transition of 
the gas (e.g. [ Elmegre en fc Parravanolll99 j : ISch ave 2004; 
iKrumholz et al.ll2009D .~ 

Outer disks have consequently been considered in- 
hospitable environments for star form ation despite the 
existence of Ha knots at large radii (Kennicuttl 119891 : 
IMartin fc Kennicuttl 120011 : iFerguson et al.lll998D . Inter- 
est in outer disks was rekindled by obse rvations with 
the G alaxy Evolution Explorer (GALEX: Ma rtin et al.l 
(2005)) which revealed that ~ 30% of disk galaxies 
have UV emitting sources beyond their optical disks 
(iThilker et aUl2007UGil de Paz et all 2005l: IThilker et al.l 
120051 : IZaritskv fc Christleinl 120071 ). In the case of M83 
and NGC 4625, UV knots have been identified as low 
mass stellar complexes and, if vi sible in the Ha, are gen - 
erally ionized by a sing le star ([Gil de Paz et al.l [2007). 
These knots are dynamically cold and rotating, indicat- 
ing t hat outer disk complexes are extensions of the inner 
d isk (IChristleiri fc Zaritskvl 120081). 

IThilker et al.l (j2007l hereafter T07) classified extended 
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UV (XUV) disks into two types. In 2/3 of XUV disks, 
the UV emission is structured, lying in filamentary, spiral 
patterns (Type I XUV disks). In the other 1/3, a large 
zone in the outer regions of the galaxy is dominated by 
UV emission (Type II XUV disks). 

The prevalence of star formation in outer disks raises 
new questions about the nature of star formation in dif- 
fuse environments. In particular, we address whether 
a threshold density to st ar formation i s consi stent with 
the observed emission. Bois sier et al.l ([20071 hereafter 
B07) selected a sample of galaxies of large angular size, 
including XUV and non-XUV disks, and studied their 
extinction corrected UV and gas profiles. While the un- 
certainties are large, B07 do not find a departure from 
the Kennicutt-Schmidt relation at low gas density. They 
suggest that the high mass end of the stellar initial mass 
function (IMF) is not sampled in every stellar cluster at 
the low star formation rates observed and Ha emitting 
stars are not formed in all star forming complexes. They 
argue that this is responsible for the difference in the 
UV and Ha profiles and the Kennicutt-Schmidt law is 
fo llowed to low surface density. 

iPflamm-Altenburg fc Kroupal ([2008) propose instead 
that the Kennicutt-Schmidt law is obeyed below the tra- 
ditional threshold density but that high mass stars capa- 
ble of emitting in the Ha simply do not form at low gas 
densities. This is supported by a deficiency of Ha emis- 
sion in dwarf gal axies that canno t be explained by IMF 
sampling effects (jLee et al.ll2009l ). iRoskar et aD (|2008at ) 
note that stars migrate radially from their formation sites 
in the disk, which can modify the stellar populations of 
outer disks. However, it is not clear that this process is 
efficient enough to rj jc-cate UV emitting stars to 2-3 D25 

Alternatively, at some sites the gas density may ex- 
ceed a star formation threshold locally, allowing stars 
to form beyond the radius where the azimuthally av- 
er aged_gas_d^nsity i s at or below a thresho l d den- 
sity ( | Kemncirt^fT9 89: Martin & Kennicutt 200 fl ISchavel 
l2QQ4HElmegreen fc Hunterll2006l : iGil de Paz et al.ll2007l ). 
Outer disk stellar complexes are often coincident with 
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local HI over-densities ([Ferguson et al.|[T998[ T07), sup- 
porting this interpretation. XUV disks are gas rich for 
their star formation rates, indicating that a high frac- 
tion of their gas is not undergoing star formation (T07). 
Studies of UV selected sources in the XUV disk of M 83 
indicate that star formation occurs only above a critical 
gas density defin ed by the Toom re criterion for dynam- 
ical instability (|Dong et al.l l2QQ8h and, where it occurs, 
follows a Kennicutt-Schmidt law. 

Thi s theory has been explored bv lElmegreen fc Hunterl 
(2006). In their model, gas clumping triggered by spi- 
ral density waves, radial variations in the interstellar 
medium (ISM) turbulence, and gas phase transitions lead 
to localized regions of active star formation in extended 
gas disks. They implement a Kennicutt-Schmidt law and 
a radially varying threshold density ba sed on the Toomre 
criterion (El megreen fc Hunterl 120061 ) . These processes 
enable star forma tion to extend into the outer disk. In 
iBush et all (|2008f ). we used smoothed particle hydrody- 
namic (SPH) simulations of an isolated disk galaxy with 
an extended gas disk of initially constant surface density 
to explore the morphology resulting from such in situ 
outer disk star formation. We showed that spiral struc- 
ture can propagate from the inner disk to the extended 
gas disk, producing regions of locally enhanced gas den- 
sity and triggering star formation in filaments similar to 
those observed in Type I XUV disks. However, this study 
was restricted to one galaxy model, limiting our ability 
to draw general conclusions. 

In t his paper, we expand on the work of IBush et al.l 
(12008ft by varying the structure of the simulated galaxy, 
in particular the gas profile and the disk to halo mass 
ratio, to explore the amount and morphology of outer 
disk star formation that result naturally without alter- 
ing traditional star formation laws. In addition, we pro- 
duce simulated images to directly compare UV and K s 
emitting populations to observations. This allows us to 
classify our models as normal, Type I or II XUV disks. 
We also compare our results to studies of the Kennicutt- 
Schmidt law, radially and locally. While we test whether 
the UV and K s emission of XUV disks can be reproduced 
using traditional star formation laws, we do not create 
Ka images and therefore do not address the discrepancy 
between Ha and UV radial profiles. In §2 we describe our 
models, in §3 we describe the evolution and uncertainties 
in one model in detail, in §4 we describe the full set of 
models and how they compare to observations and finally 
in §5 and §6 we discuss our results. Our models succeed 
in producing the UV morphologies of Type I XUV disks. 
However, they do not reproduce the UV bright outskirts 
of Type II XUV disks. 

2. METHOD 

2.1. SPH Simulations 

We use the code Gadget2 (|Springell 120051 ). based on 
the fully conservative formu lation of SPH developed by 
ISpringel fc Hernquistl ()2QQ2h to run simulations of galax- 
ies with extended gas disks. We incorporate a sub- 
resolution multi-p hase model of the ISM and s tar forma- 
tion according to Spr ingel fc Hernquistl (|2003 | . hereafter 
SH03) , whi ch includes r adiati ve cooling as in lKatz et al.l 
(1996) and I Dave et al.l ([1999). Star formation from the 
cold phase gas (total cold gas, our simulation does not 
treat atomic and molecular gas separately) follows a 



Schmidt volume density law Psfr oc p^ as with TV = 1.5 
and normalized to a star formation timescale (£*). 

We run simulations with two different star formation 
timescales to explore how t* affects the XUV emission. 
One set of simulations adopts the values proposed in 
SH03: U = 2.1 Gyr, A = 1000, and T SN = 10 8 , 
where Ao and Ts n are parameters that regulate the mass 
fraction in the cold phase of the ISM. This normalizes 
the star fo r mation rate to the best-fit relation found in 
IKennicuttl (|l998bl ). The second set of simulations are 
run with parameters which normalize the star formation 
rate roughly to that of the Milky Way (MW), since these 
are MW like models: t* = 6.25 Gyr, A = 3000, and 
Tsn — 3 x 10 8 . This results in lower levels of star for- 
mation for a given gas density. The local star formation 
volume density cutoff corresponds to ~ 0.004 M /pc 3 
for both sets of parameters. In all cases, s tar forma- 
tion rat es lie well within the scatter shown in IKennicuttl 
(1998b) and the surface density threshold fal ls between 
~ 3-5 Mfp/pc 2 , consistent with observations ([Kennicuttl 
119891 : iMartin fc Kennicuttll200l[ ). The effective equation 
of state for the star-forming gas is determined from these 
parameters through the pressurizing action of supernova 
feedback, allowing us t o stably evolve disks with even 
very large gas fractions ([Springel fc Hernqu ist 2005). 

2.2. SPH Initial Conditions 

Our galax y models are constru cted using the procedure 
outlined in ISpringel et al.l (j2QQ5h . We construct analogs 
to the Milky Way with a total mas s of M200 ~ 10 12 
M composed of a iHernquistl (|l990l ) dark matter halo 
with concentration c = 9, an exponential stellar disk 
with scale length 3.75 kpc (we adopt the Hubble pa- 
rameter h = 0.72 throughout this paper), and an ex- 
ponential gas disk of the same scale length. We then add 
an extended gas disk either of appr o ximat ely constant 
density or following an r _1 Mestel (1963) radial pro- 
file that starts at a "break radius" of 3.5 scale lengths 
(~ 13 kpc) and extends to 8.7 scale lengths (~ 32 
kpc). We also create a control model with no exten- 
sion (pure exponential profile). We then normalize with 
each of these outer disk profiles (pure exponential, fiat, 
and Mestel) to have a gas surface density of 3, 5 and 
10 M pc _2 at the break radius, spanning the range of 
outer disk profiles seen in the THINGS (jWalter et alJ 
200^) and lBroeils"fc van Woerdenl (jl994f ) surveys. Ta- 
ble [T] summarizes these properties and lists a run-name 
for each model. We do not simulate a galaxy with 
a high mass constant density extension because nearly 
constant density disks are not found at those densi- 
ties. The gas content and distribution agree roughly 
with the Mm/Mtot an d log Mm / log Dm trends seen in 
iBroeils fc Rheel (|l997f ). T he fiat-low model is very simi- 
lar to the model shown in lBush et al.l (|2008l ). 

The ratio of the mass of the inner gas disk to the mass 
of the inner stellar disk is fixed at ~ 1:5, meaning the 
overall stellar disk mass is increased when the gas surface 
density is increased. The halo mass is also increased ac- 
cordingly to keep the halo to disk mass ratio at constant 
at 0.041. Varying the halo to disk mass ratio changes the 
spi ral structu r e amp litude, which we do in § 15.11 A stel- 
lar [Hern^uisa (|l990l ) bulge with a scale length of 0.27 kpc 
is added with 0.008 of the total mass of the simulation, 
making it ~ 25% the mass of the stellar disk. Particle 
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TABLE 1 

Gas distributions for models used in this study. Every 
model is initialized with an exponential inner disk profile 
out to 13 kpc and the outer disk profile listed here to 32 
kpc. The mass of each component of the galaxy is scaled 

to give the surface density listed here at the break 
radius and the number of particles in each component is 

scaled to keep the mass per particl e co nstant in each 
component. Details are described in § 12.21 All models are 
evolved in isolation. 



numbers are varied for each model to keep the mass per 
particle constant. The mass per particle in the stellar 
and gas disks is ~ 1O 5 M and in the halo is ~ 1O T M . 
The Gadget softening length of the baryonic particles 
and halo particles are 0.14 kpc and 0.28 kpc, respectively. 
The models each have a 0.01 M200 mass black hole in its 
center. The black hole is represented by a sink particle 
that can accrete mass at th e Bondi rate ([Springel et al.l 
l2QQ5t iDi Matteo et al.ll2005h . but as we are interested in 
the outer disk its presence does not influence our results. 

2.3. Radiative Transfer 

To create simulated FUV and K s -band images, we 
employ the 3D adaptive grid polychr omatic Monte 
Carlo radiative transfer code Sunrise 5 (|Jonssonl [20061 : 
Uonsson et aD l2009h . SUNRISE assigns a spectrum to 
each particle in our simulations that would be a source 
of radiation (stellar and black hole particles) and then 
propagates photon "packets" from these sources through 
the dusty ISM, represented by the gas particles, using 
a Monte Carlo approach. The black hole particle emits 
an empirical template AGN spectrum der ived from ob- 
servations of quasars (Ho pkins et al.l 12007) but since we 
are examining the outer disk, and this is a small black 
hole with almost negligible accretion, this has negligible 
impact on our results. 

Each stellar particle over 10 Myr old is given a spec- 
trum corresp onding to its age a nd metallicity from 
Starburst 99 ([Leitherer et al.l fl999[ ) populati on synthe- 
sis m odels. We assume a Kroupa IMF ([Kroupa et al.l 
Il993f ). Star particles that existed at the beginning of 
the simulation are assigned an age from a distribution 
that is determined by a star formation history chosen to 
match observations of local disks (see lRocha et al.l {2008) 
for more details). The disk was assigned an exponen- 
tially increasing star formation rate with a timescale of 
106 Gyr (which makes the star formation rate effectively 
constant) with the oldest stars being 13.9 Gyr old. The 
bulge was assumed to form in an instantaneous burst 
with an age of 13.9 Gyr. Since these parameters are not 

5 Sunrise is a free, publicly available software that 
can be applied to a ny hydrodynamic galaxy simulation, 
http: / / sunrise.familjenjonsson.org 




Fig. 1. — Metallicity as a function of radius in our galaxy models 
at 0.14 Gyr (initial snapshot). The solid line describes the stars 
and the dashed line follows the gas metallicity into the outer disk. 



given any radial gradient, and the number of stars in the 
outer disk is low, they only slightly affect our results. 

Star clusters with ages less than 10 Myr are assumed 
to be located in their nascent birth-clouds of molecu- 
lar gas and are given a modified spectrum which ac- 
counts for the effects of H II and photo-dissociation re- 
gions (PDRs). The evolution of HII regions and PDRs 
are described anal ytically and the ph oto-ionization code 
MAPPINGS III ([Groves et al.ll2004D is used to calculate 
the propagation of the s ource spectrum through this neb- 
ula {G roves eTaIll200cl . The HII region absorbs effec- 
tively all ionizing radiation and is the source of hydrogen 
recombination lines. A fraction (/pdr=-3) of the solid 
angle into which the star particle is emitting is assumed 
to be obscured by the PDR, which absorbs a significant 
fraction of the UV emission of the stellar particle and 
reprocesses it into the FIR. The remaining solid angles 
are not influenced by the PDR, and the final result is 
a linear combination of the PDR of obscured and unob- 
scured regions weighted by /pdr- Owing to the low star 
formation rates in the outer disk, very few particles in 
the outer disk are less than 10 Myr old. Thus emission 
from stars unaffected by PDRs dominates the FUV and 
the results here are very insensitive to /pdr- 

The source spectra, as described above, are then prop- 
agated as photon 'packets' through the ISM represented 
by gas particles. Gas and star particles are initialized 
with an exponentially decreasing metallicity distribution 
with a central value of 1.6 Z and a value of 0.3 Z at 
the break radius to a pproximately match the outer disk 
metallicities of M83 dBresolin et al.l l2QQ9h . The distri- 
bution is shown in Figure Q] and additional metals are 
added as the simulation progresses as described in SH03. 
We a ssume a dust to metals m ass ratio of 0.4 dDwekl 
I1998T) and use the dust models of Weingartner fc Draind 
(|200lD . updated according to iDraine fc Lil (120071). Fo r 
details on the photon propagation, see lJonssonP j2006). 
The final spectra are comprised of source photons that 
escape absorption and scattering, supplemented by pho- 
tons re-emitted in the infrared from the dust and photons 
scattered into the line of sight. Finally, the viewing angle 
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and bandpass are considered to calculate broadband im- 
ages. Since we examine only FUV and K s -band images in 
this paper, dust emission is unimportant. We only exam- 
ine face-on images in this work. Highly inclined models 
would have more obscuration along the line of sight and 
show less UV emission. 

3. A CASE STUDY 

3.1. Evolution 

In lBush et~a l. (2008) we demonstrated that spiral over- 
densities in the inner disk of a galaxy can propagate to 
the edge of an extended constant density gas disk. In the 
outer disk of that model, peaks in the spiral waves above 
the density threshold form stars while troughs below the 
density threshold do not. This yields filamentary star 
formation in the outer disk, as seen in Type I XUV disks 
(T07). A study of the model's observable properties us- 
ing radiative transfer was beyond the scope of that work. 
However, the model presented there is very similar to our 
flat-low model. To provide a context for the analysis of 
the entire sample of simulations, we first study the re- 
sults of applying radiative transfer to the flat-low model 
and the time evolution of this model's properties. 

Four snapshots from the 4 Gyr evolution of flat-low are 
shown in Figure [2] The first three rows of Figure [2] map 
stellar surface density, gas surface density, and instanta- 
neous star formation rate surface density (top-bottom). 
The star formation rate, and therefore also the UV emis- 
sion, is determined by the gas distribution. Early in the 
simulation, ring-like perturbations in gas density prop- 
agate outwards from the center. These relax with time 
into steady spiral features caused by potential fluctua- 
tions in the halo (He rnqui"stl Il993) . For a discussion of 
how these spiral features compare to observed spiral fea- 
tures in galaxies, see § 15.11 This results in rings of star 
formation in the outer disk early on, relaxing to spiral 
features around 2 Gyr. Late in the simulation, the inner 
regions deplete as gas is converted into stars, leading to 
inner holes in the gas distribution and, eventually, also 
in the distribution of star formation. The star forma- 
tion rates in the outer disk are low enough that they do 
not significantly deplete the outer disk (see § 15. 4ft . so the 
morphology of star formation in the outer disk changes 
very little from 2 Gyr to the end of the simulation. 

To classify disks according to their outer disk star 
formation, T07 analyze the FUV and 2MASS K s -band 
images of 189 galaxies in the G ALEX Atlas of Nearby 
Galaxies (iGil de Paz et al.ll2007f ) and define two types of 
XUV disks. FUV emission traces star formation over ~ 
200 Myr, not instantaneous star formation, so in order 
to accurately compare our simulations to T07 we cre- 
ated FUV and K s -band images with Sunrise. FUV and 
K s -band images are shown in the fourth and fifth rows 
of Figure O The FUV emission partially traces the in- 
stantaneous star formation rate, but also stretches the 
features out in time as older stars continue to emit in 
the UV. The K s -band traces the complete stellar pop- 
ulation, resulting in images very similar to the stellar 
surface density maps. 

T07 define a Type I XUV disk to be a galaxy that has 
at least two complexes of structured FUV emission be- 
yond an FUV contour of 27.25 AB magnitudes arcsec -2 
at 1 kpc resolution that are unaccompanied by emission 
in the K s or DSS red bands. To analyze our simulations 



when the gas distribution is comparable to observed spi- 
ral galaxies, we choose snapshots after the spiral struc- 
ture has become steady, but before the inner gas is de- 
pleted, near the 2.285 Gyr snapshot shown. At this time, 
knots of FUV emission trace spiral arms in the outer disk. 
The features brighter than 27.25 AB mags arcsec -2 are 
highlighted with a thin black contour in Figure [2] The 
surface brightness of the outer disk in the K s is very low. 
The approximate surface brightness limit of 2MASS is 
20 magnitude s arcsec -2 , shown as the white dashed line 
in Figure l2l ( Jarrett et al.ll2003[ ). Most of the emission in 
the outer disk is below the detection limit for 2MASS. 
The outer disk shows more than one knot of UV emis- 
sion without accompanying K s -band emission, this model 
most likely fits the criterion for Type I XUV disks. We 
do not simulate the DSS red band, but expect the results 
to be similar. 

T07 define Type II XUV disks based on the proper- 
ties of a low surface brightness (LSB) zone in the outer 
disks of galaxies. They define the LSB zone as emis- 
sion beyond 80% of the K s -band flux, and interior to the 
innermost 27.25 FUV AB magnitudes arcsec -2 contour. 
The 80% K s -band flux contour for our model (80 % of the 
flux above the 20 magnitudes arcsec -2 surface brightness 
limit) is shown in Figure [2] in yellow. The innermost con- 
tinuous 27.25 FUV mag arcsec -2 contour is highlighted 
in green on the same figure. The LSB zone of our sim- 
ulations is then the area between the yellow and green 
contours. T07 use the color and area of the LSB zone to 
define Type II XUV disks as disks whose outskirts are 
very bright in the FUV, but low surface brightness in 
the K s -band. Their empirical definition of these objects 
is those whose LSB zones have an FUV— K s color of less 
than 4, and whose LSB zones are at least 7 times the area 
of the inner disk, defined by the 80% K s -band contour. 

For this model, at 2.285 Gyr, the LSB zone color is 
3.6, but the LSB zone size ratio is 2.1. Consequently, 
flat-low does not qualify as a Type II disk at this time. 
In fact, these values are similar to those sees in Type I 
XUV disks and normal disk galaxies (see § !4.ip . Because 
star formation in the outer disk is very filamentary, the 
continuous FUV contour, even at 1 kpc resolution, ac- 
tually traces the edge of the inner disk. Therefore, the 
LSB zone is actually quantifying the properties of the 
edge of the inner disk, not the properties of the extended 
outer disk, so we would expect its properties to agree 
with those of normal galaxies. 

Since the LSB zone does not measure extended fila- 
mentary emission, we introduce a covering fraction of 
FUV emission to quantify the amount of outer disk emis- 
sion. We define the covering fraction as the fraction of 
pixels beyond the 'FUV inner disk' that are brighter than 
27.25 FUV AB magnitudes arcsec -2 when the image is 
smoothed to 1 kpc resolution. We define the FUV inner 
disk to end where the average FUV surface brightness 
is 27.25 FUV AB magnitudes arcsec -2 and calculate the 
covering fraction between this radius and twice this ra- 
dius. The values are listed below the FUV images. These 
will be discussed further in § [U 

To explore how the evolution of the disk affects these 
classifications, we examine the properties of the outer 
disk with time. We take the outer disk to start at a ra- 
dius where the star formation rate surface density falls 
below 3 x 10 -4 M kpc -2 yr -1 , the 27.25 FUV contour 
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Fig. 2. — Snapshots of the evolution of the flat-low model. Top-bottom: stellar surface density, gas surface density, star formation rate 
surface density, FUV surface brightness and K s surface brightness. The black contour is the 27.25 FUV AB magnitudes arcsec -2 contour 
and the yellow contour contains 80% of the K s -band flux. The green contour highlights the innermost 27.25 FUV AB magnitudes arcsec - 2 
contour. The area between the yellow and green contours is defined as the LSB zone in T07. The CF label is the covering fraction of outer 
disk FUV pixels above 27.25 FUV AB magnitudes arcsec -2 . The white dotted line is at 20 K s Vega magnitudes arcsec -2 . 



in star formation rate using the calibration in Kennicutt 
(1998a), when the simulation is initialized. The mass of 
stars initially in the simulation, gas mass, star forma- 
tion rate and mass of stars formed in the simulation are 
summed beyond this radius every 0.14 Gyr and presented 
as a function of time in Figure [3l For the first 1.5 Gyr, 
as the simulation relaxes, strong oscillations are seen in 
these quantities as density perturbations shift slightly in- 
side or outside the boundary between the inner and outer 
disk. After these relax, most quantities evolve smoothly 
with time. The mass of old stars in the outer disk in- 
creases slowly over time. This most likely owes to stellar 
migration, a phe nomenon caused by resonant scattering 
from spiral arms (|Roskar et al.ll2QQ8aj ) . The gas mass de- 



creases slowly with time as it is converted into stars. The 
mass of stars formed increases correspondingly. The star 
formation rate is high initially owing to the strong pertur- 
bations in gas density as the disk relaxes, but then slowly 
decreases as gas depletes. As the stellar mass of the disk 
builds and the star formation rate drops, the color of the 
LSB zone reddens and the size shrinks. Since Type II 
XUV disks have large, blue LSB zones, this means that 
this disk evolves away from being a Type II XUV disk 
over time. This hints that recent gas accretion, creating 
in effect an un-evolved outer disk, may be important in 
producing Type II XUV disks. The covering fraction of 
FUV emission in the outer disk is also high initially but 
settles at around 0.02 after 2 Gyr. It is noisier than many 
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Fig. 4. — Average metallicity of the flat-low model as a function 
of radius for the stars formed in the simulation (solid lines) and 
gas (dashed lines) as a function of time. Times are represented by 
colors: 0.7 Gyr - purple, 1.4 Gyr - blue, 2.1 Gyr - turquoise, 2.8 
Gyr - green, 3.5 Gyr - yellow and 4.2 Gyr - red. 

The results for the double resolution run are over-plotted 
in Figure [3] in orange. The observable quantities do not 
change more than the scatter in the original simulation. 
Although the simulation with ten times the resolution 
simulation was too large to be processed by Sunrise, we 
have over-plotted the stellar mass, gas mass and star for- 
mation rates of the outer disk in Figure [3] in green. Note 
that at this resolution, the star formation rate increases. 
The simulation is now tracing smaller, but higher den- 
sity peaks in the gas density distribution. We note that 
this means our results are lower limits on the star forma- 
tion rate, and most likely the covering fraction, of outer 
disks. However, the qualitative results described in § [4] 
are unaffected. 



Fig. 3. — Time evo lution of outer disk quantities for the flat-low 
model. See the § 13. II for a description. 

of the other quantities, but evolves less with time. 

To confirm that these simulations are indeed relaxed 
after « 2 Gyr, we created new initial conditions for the 
flat-low model by pre-running the simulation without any 
star formation or gas cooling for over 4 Gyr. We then 
turned star formation and cooling on and the simulation 
was analyzed identically to the original. After 2 Gyr, the 
results are consistent with the original run. 

3.2. Dependence on Disk Resolution 

The three observational quantities (covering fraction, 
LSB zone color and LSB zone size) are evaluated over an 
area that depends on the emission properties of the disk. 
If the area covered by star forming particles changed sig- 
nificantly with resolution, these quantities may be reso- 
lution dependent. The covering fraction, with its depen- 
dence on the knots of star formation in the outer disk, 
is particularly susceptible. To determine whether these 
quantities depend on resolution, we ran simulations with 
twice and ten times the number of particles in the gas, 
stars and bulge of the galaxy. In these simulations, more 
gas particles of lower mass, form stars at each time step. 



3.3. Dependence on Metallicity 

The metallicity of the stars in our simulation affects 
the template SED chosen for them, and the metallicity 
of the gas affects the amount of obscuration in the ra- 
diative transfer calculations. It is therefore important 
to consider how the metallicity distribution we select af- 
fects our results. The metallicity used by Sunrise comes 
from two sources. The first is the distribution of metals 
present when the simulation is initialized. Our initial 
distribution is shown in Figure [TJ The second is metal 
enrichment owing to supernova feedback. As stars form 
in the simulation, gas particles are treated as a closed 
box system, and metals are instantaneously returned to 
the gas. The new stellar particle is given the metallicity 
of the gas particle from which it formed. Over time, the 
metallicity of the gas and stars increases. This evolution 
is shown for the flat-low simulation in Figure HI 

Very little is known about the metallicity of 
outer disks. Observations have found gas metallic- 
it ies anywhere from 1 — .1 Z ( Bresol in et al.l 2009; 
iGil de Paz et al.l 120071 ). At 2.2 Gyr, the metallicity of 
the new stars (the dominant source of UV emission) at 
the break radius in our simulation is ~ 0.3 Z from su- 
pernova feedback and ~ 0.3 Z from the initial condi- 
tions, which gives ~ 0.6 Z overall. This is reasonable, 
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Time (Gyr) 

Fig. 5. — Time evolution of outer disk quantities for the flat-low 
model. Over-plotted in orange are the results for a simulation with 
no initial metals. See the § 13. II for a description. 

but could also be too high. To test how much our results 
depend on metallicity, we re-ran the Sunrise radiative 
transfer on flat-low with no initial metal content. In 
this case, the entire metal contribution comes from su- 
pernova feedback. Clearly this is a limiting case, since 
it also assumes that the initial simulation stars have no 
metal content. The results are shown in Figure [5] The 
black line is the flat-low simulation and the orange line 
shows the results of the radiative transfer with no ini- 
tial metal content. As we would expect, the LSB zone 
color is noticeably bluer throughout the simulation, ow- 
ing to the lower metallicity of the FUV emitting stars. 
The LSB zone size increases owing to an increase in FUV 
emission and a slight decrease in K s -band emission, both 
resulting from the lower metallicity of the initial and new 
stars. Interestingly, the covering fraction changes very 
little. This indicates that the change in metallicity does 
not increase the brightness of knots enough for a larger 
number to exceed the 27.25 FUV contour. 

While these effects are not large enough to change the 
conclusions we discuss, there is some uncertainty about 
the exact color and sizes of the LSB zone owing to the 
lack of constraints on the metallicities of outer disks. 

4. SET OF MODELS AND COMPARISON TO 
OBSERVATIONS 

After the detailed examination of the flat-low model, 
we now study and classify the other models in Table HJ 
Spiral structure propagates from the inner to the outer 
disk in all of our simulations. The amount of outer disk 
star formation is then controlled by the amount of gas in 
the outer disk that is above the threshold at any given 
time, which depends on the density perturbations owing 
to the spiral structure and the radial distribution of the 
gas. Immediately below, we focus on the radial gas distri- 
bution, but we consider variations in the spiral structure 
in §0 

Maps of stellar surface density, gas surface density, and 
instantaneous star formation rate, along with gas, stellar, 
and star formation rate density profiles at 2.2 Gyr for 
each simulation are provided in Appendix A. The time 
evolution of each model is similar to that seen for the flat- 




CF=0.006 ±0.006 CF=0.002±0.001 CF=0.003±0.002 




CF=0.009 ±0.006 CF=0.008 ±0.006 CF=0.006±0.005 




CF=0.020 ±0.007 CF=0.129±0.049 



Fig. 6. — 2.2 Gyr FUV broadband image from each of our sim- 
ulations. Top (L-R): exp-low, exp-mid, exp-high. Middle (L-R): 
Mestel-low, Mestel-mid, Mestel-high. Bottom (L-R): flat-low, flat- 
mid. Contours are described in Figure [2] 




Fig. 7. — K s -band images for the same simulations as in Figure[6] 
The white dotted line is at 20 K s Vega magnitudes arcsec - . 

low model. In the following sections, we compare our 
models to several observations star formation in outer 
disks. 

4.1. Comparison to Thilker et at. 2007 

One snapshot from each simulation described in § 12.21 
is shown in Figure El The thin black contour in these im- 
ages is at 27.25 AB magnitudes arcsec -2 , and the green 
line highlights the innermost continuous contour, as in 
Figure [2] Clumps beyond this green line contribute to 
the Type I XUV disk classification. All disks show a few 
clumps beyond this contour, but how fast the gas distri- 
bution falls off with radius controls the radial extent of 
the FUV complexes. Disks with slowly varying gas dis- 
tributions have clumps of FUV emission at larger radii. 
The flat-low, flat-mid, Mestel-mid, and Mestel-high mod- 
els all show significant numbers of FUV knots in filamen- 
tary structures extending off the disk. None of the purely 
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CF=0.002 ±0.001 CF=0.002 ±0.002 CF=0.004 ±0.003 




CF=0.009 ±0.004 CF=0.010±0.005 CF=0.007±0.006 




CF=0.009 ±0.002 CF=0.089 ±0.039 



Fig. 8. — Same as Figure [6] for simulations with t* = 6.25. 

exponential gas disks show more than a small 'fringe' of 
XUV emission because of their rapid decrease in surface 
density with radius. The corresponding K s -band images 
are shown in Figure [71 The surface brightness limits seen 
in Figure [7] (white dashed contour) are at smaller radius 
than most of the FU V knots seen in Figure [6j so the flat- 
low, flat-mid, Mestel-mid, and Mestel-high models would 
most likely be classified as Type I XUV disks. 

To determine how robust these results are to the star 
formation rate prescription we increased the star forma- 
tion timescale of the gas by a factor of ~ 3. Shown in 
Figure [8j when t* is increased, the overall amount of star 
formation is reduced, but its distribution is largely un- 
changed. However, the flat-low model and Mestel-mid 
model simulation may not be classified as XUV disks, 
depending on the snapshot analyzed, because they have 
very few knots of emission brighter than the 27.25 FUV 
contour. 

Based on these simulations we expect higher surface 
density gas disks whose outer gas profiles do not fall 
off exponentially to be classified as Type I XUV disks. 
Studies of H I distributions indicate that non-exponential 
gas profiles are comm on in the outer parts of gala xies 
(|Broeils fc Rheel 1199ft iBroeils fc van Woerdenl 119941 ). It 
should be noted that we only tested exponential disks 
with the same scale length as the stellar disk, in this 
case, 3.75 kpc. As the scale length of the gas disk in- 
creases and falls off more slowly it will be more likely to 
be classified as an XUV disk. 

To determine whether any of the models would be clas- 
sified as a Type II disk, we evaluate the color and area 
of the low surface brightness (LSB) zone defined by T07. 
As described in § 13.11 T07 define the LSB zone as emis- 
sion beyond 80% of the K s -band flux, and interior to the 
innermost 27.25 FUV AB magnitudes arcsec -2 contour. 
The 80% K s -band flux contour for the grid of models is 
shown in yellow in Figure [6l The LSB zone of our grid 
of simulations is the area between the yellow and green 
contours in Figure [6l T07 define Type II XUV disks as 
objects whose LSB zone has an FUV— K s color of less 
than 4, and the LSB zone is at least 7 times the area of 
the inner disk, defined by the 80% K s -band contour. The 
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Fig. 9. — LSB zone diagnostic diagram to compare to T07, with 
the ratio of the area of the LSB zone to the area containing 80% 
of the K band flux vs. the K s - FUV color of the LSB zone. For a 
description of the LSB zone see the text. T07 data are represented 
by black points for non-XUV disks, red squares for Type I XUV 
disks, and blue diamonds for Type II XUV disks. Our models are 
shown as crosses and models we identified as Type I XUV disks are 
circled. The most high mass models are in blue, middle in green 
and low mass in orange. Points are the average values over four 
snapshots and error bars are the standard deviation of these values. 
Our models are consistent with being normal or Type I XUV disks. 
The inset shows the detail of the region where our models lie. 

color vs. area of this zone for our sample are shown in 
Figure [9] with the T07 sample. 

None of our simulations meet the Type II criterion. Al- 
though they are blue enough in FUV-K S , the LSB zone is 
too small. These characteristics are consistent with nor- 
mal and Type I disks. Note that our simulations have 
very small scatter in this plot. As with the discussion 
of the flat-low model, the outer disk star formation in 
our simulations covers such a low fraction of the disk 
that the innermost continuous 27.25 FUV AB magni- 
tudes arcsec -2 contour tends to fall at the edge of the 
inner disk, not the outer extended gas disk, so the LSB 
zone is primarily probing a part of the disk where the 
initial structure of each model is very similar. By vary- 
ing the structure of the inner disk, particularly increas- 
ing the scale lengths of the models, we could change the 
characteristics of the LSB zone. However, it seems un- 
likely that we could alter the structure of the inner disk 
enough that it would yield Type II disks, while staying 
consistent with the observed structure of galaxies. 

To quantify the outer disk UV emission we calculate 
covering fractions as described in § 13.11 Since the cov- 
ering fraction varies significantly between snapshots, we 
considered 4 snapshots spanning an evolutionary time 
of 2.1 Gyr to 2.5 Gyr and report their average value in 
Figure El The error reported is the standard deviation 
of these values. Overall, the covering fraction seems to 
quantify the amount of FUV emission well, but, espe- 
cially when limited to such a small sample, the covering 
fraction does not differentiate strongly between galaxies 
with different gas distributions. On average, the Mestel 
models have a higher covering fraction than the exponen- 
tial models, and likewise for the constant density models 
vs. the Mestel models. But the covering fraction is not 
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Fig. 10. — Left: SFR vs. HI mass for our simulations at 2.2 Gyr (symbols as in Figure [9]) compared to T07. Our models are consistent 
with normal and Type I XUV disks. Right: SFR/H I vs. Lf~ for our simulations at 2.2 Gyr (symbols as in Figure[9]). Our models classified 
as Type I XUV disks appear to lie in a region mostly populated by Type I XUV disks. 



significantly higher for the models classified as Type I 
XUV disks than for the others. This at least partially 
owes to resolution. The contours are smoothed to 1 kpc 
in all models, the area of the low mass disks is smaller, 
increasing the size of the UV knots in the outer disk rel- 
ative to the area the covering fraction is calculated over. 
This biases the covering fractions to be larger for the 
lower mass disks. 

To further constrain our models, we compare their 
properties to other observational characteristics of XUV 
disks. T07 show that XUV disks are gas rich for their 
total star formation rate. In Figure [10] we plot the ap- 
proximate H I gas masses of our disks against their total 
instantaneous star formation rates, along with the T07 
sample. The H I gas mass for our simulations was ap- 
proximated by a ssuming; that 50% of all the star forming 
gas is molecular ([Keres et al. 2003) and that 70% of what 
remains is HI (to account for He). Our models all lie in 
the high gas mass region of the T07 data, unsurprising 
since we have chosen a massive galaxy as our fiducial 
model. Additionally, in star forming regions, the molec- 
ular gas percentage may in fact be much higher, lowering 
the HI masses. However, since they are high mass disks, 
all of our galaxies lie in areas where there are both Type 
I XUV disks and normal disks. So while our results are 
consistent with T07, this plot does not provide a good 
diagnostic for determining whether our models match ob- 
served XUV disks in terms of gas richness. 

In Figure [10] we also plot the SFR/H I ratio vs. the 
K s -band luminosity (calculated inside the 20 magnitudes 
arcsec -2 contour) of our galaxies, and compare these to 
T07. This plots gas richness vs. Lk, a proxy for stellar 
mass and show whether the disks are gas rich for their 
star formation rate and stellar mass. Our modeled Type 
I XUV disks fall in an area that is dominated by Type I 
XUV disks, rather than normal galaxies. This is promis- 
ing, indicating that the structure of our model galaxies 
is similar to Type I XUV disks. 

In each of these plots it is evident that our models 
with constant density extended gas disks tend to have 



slightly larger HI masses than most observed galaxies. 
This is unsurprising; in reality we would expect to find 
few galaxies whose H I disks do not fall off at all with ra- 
dius. However, they are an interesting limiting case. Our 
models with r _1 extended disks seem to be quite consis- 
tent with the HI masses of massive XUV disks. Overall, 
the structure of our models that are classified as Type I 
XUV disks are consistent with XUV disks observed by 
T07. 

It is difficult to determine the fraction of XUV disks 
that owe to isolated in situ star formation if we do not 
know the fraction of XUV disks that have disk structures 
similar to that of our Type I XUV disk models. De- 
tailed studies of the H I profiles of XUV disks are needed 
to determine whether we are describing them fully or if 
there are additional conditions needed to create all Type 
I XUV disks. 

4.2. Comparison to Boissier et al. 2007 

B07 studied the extinction corrected UV radial profiles 
of 43 galaxies which were selected for large angular size. 
In contrast to studies of Ha profiles of galaxies, which 
find a stron g truncation at 3-5 M^pc~ 2 in total gas sur - 
face density (|Kennicuttlll989l : lMartin fc Kennicuttll200ll ). 
B07 find that star formation rate profiles derived from 
UV profiles do not show a break. In Figure [TT] we show 
the instantaneous star formation rate of our simulations 
vs. the gas surface density. UV derived star formation 
rates give the star formation rate averaged over 200 Myrs, 
but since the star formation rates of our simulations are 
mostly constant with time after relaxing initially, the 
trends seen should be identical. 

While at low gas densities, simulations can differ con- 
siderably in the details, all simulations, including Type I 
XUV disks, show strong truncations as a result of the vol- 
ume density threshold to star formation and azimuthal 
averaging. The levels of star formation at low gas den- 
sities indicated by B07 are simply not reflected in our 
simulations. However, there are large uncertainties in 
the low star formation rate measurements and, when a 
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Fig. 11. — Plots of star formation rate vs. gas surface density radially binned for all our models: high surface density in blue, middle 
surface density in green and low surface density in orange. Outer disk profiles are separated by symbols: exponential - plus signs, Mestel 
- triangles and constant density - diamonds. On the left, quantities are averaged in radial bins 1 kpc wide and on the right 5 kpc wide. 
While this can change the detaile d nature of the profile, it does not eliminate the truncation in the profile. The dashed line is the empirical 
law denned by Kennicutt (1998b). Data from B07 are plotted as black points. 



gas profile is slowly varying, averaging can play a sig- 
nificant role in smoothing the star formation profile. To 
demonstrate this we plot the gas density vs. star for- 
mation rate density radial profile with radial bins of two 
different sizes, 1 kpc and 5 kpc. While with larger bins 
the profiles do flatten, they still show significant trun- 
cation. If the B07 results hold true with more detailed 
studies of individual objects, a change in the star for- 
mation prescription should be needed at low gas surface 
densities for at least some galaxies. 

If the resolution of our simulation was increased, it 
is possible it could trace smaller overdensities and in- 
crease the amount of star formation at low average 
gas densities. However, when the azimuthally averaged 
Kennicutt- Schmidt law is plotted for the fiat-low simula- 
tion at ten times the disk resolution, a truncation is still 
very evident. 

4.3. Comparison to Dong et al. 2008 

iDong et al.l (|2QQ8f ) studied the characteristics of UV se- 
lected sources in the outer disks of a prototypical Type 
I XUV disk, M83, using circular apertures with a diam- 
eter of 440 - 600 pc. They found that most of the star 
forming complexes are coincident with H I gas above the 
Toomre dynamical instability threshold. In addition, us- 
ing UV and 8 jiva derived star formation rates and plot- 
ting these against HI gas density, they found general 
agreement with the Kennicutt- Sc hmidt law derived fo r 
regions in the inner disk of M 51 ([Kennicutt et al. 2007). 

In order to compare to this type of analysis, we divide 
the star formation rate and gas surface density maps of 
each simulation into 1 kpc x 1 kpc squares and plot the 
star formation rate vs. gas surface density fo r squares in 
the ou ter disk in Figure [12j along with the IDong et al.l 
(|2008l ) data. The SPH smoothing length is typically 1 
kpc for our outer disk particles, therefore these are the 
smallest apertures we can use and t hey should be com - 
pared to the smaller aperture data of IDong et al.l {2008) 
with caution. With the exception of a few points (and 
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Fig. 12. — Plots of star formation rate vs. gas surface density in 
lxl kpc apertures over the outer disk. Green points indicate star 
formation at that rate and surface density beyond approximately 
1.5 optical radii as well as between the optical radii and 1.5 optical 
radii. Red points indicate star formation and are interior to 1.5 
optical radii. Models are the same as in Figure [6] The over-plotted 
black crosses are measured SFRs and H I gas densities (lower limits 
on th e total gas density) from the outskirts of M83 (|Dong et al.l 
[20081 V 

IDong et al.l (|2008l ) caution that the uncertainties on the 
8 fim portion of the star formation rates are large), we 
create star formation rates as high in the outer disk as 
those measured in M 83. Some gas surface densities in the 
IDong et "all ([2008) data are considerably lower than our 
simulations, but they only include HI, so these should be 
considered lower limits on the gas surface density. If we 
calculate the HI surface density for our simulations as- 
suming 50% of all star forming gas is molecular (clearly 
not accurate on large s cales in out e r disk s but in the 
small apertures used by IDong et all ((2008) this may be 
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a good approximation) we find that the central concen- 
tration of the M 83 data agrees well with our simulation 
results. In addition, some o f our star formatio n rates are 
lower than those see in in iDong et all (|2008ft . but this 
is expected because they lay apertures only on sources, 
while we simply bin our outer disk, so some of our bins 
will have very low star formation rates. 

As the star formation prescription is the same in each 
model, the Kennicutt- Schmidt law averaged on local 
scales change very little as the disk structure is changed 
between models. The data agree similarly with all the 
models we considered. The main difference between the 
models with and without Type I XUV emission is the ex- 
tent of the star formation. To highlight this, we have col- 
ored portions of the diagram where star formation exists 
beyond a distance of approximately 1.5 optical radii in 
green. It is clear that the exponential disks do not show 
any star formation at large radii and, at large radii, only 
Mestel-mid, fiat-low and flat-mid agree with the data. 

5. DISCUSSION 

5.1. Spiral Structure 

To determine whether we are accurately reproducing 
the star formation observed in galaxies, we need to deter- 
mine whether the spiral perturbations in our simulations 
are consistent with those in actual galaxies. Though we 
know that spiral structure in o ur simulations is caused 
by numerical noise in the halo ([Hernquistl 1993), what 
causes spiral structure in nature is unknown. However, 
for our purposes we only need to accurately recreate the 
amount of gas above the star formation threshold and 
therefore the amplitude of the spiral perturbations in the 
outer disk. 

Studies of the spiral structure in stellar disks indicate 
that the two armed mode of spiral structure (m = 2) 
dominates in most galaxies. We plot the relative ampli- 
tude of the stellar m = 2 mode for the 2.2 Gyr snapshot 
of each simulation on the left in Figure [13] The very high 
amplitudes in inner regions owe to central bars. These 
are generally consi stent with the range of observed m = 2 
mode amplitudes (Rix fc Zarit skv 1995: lElmegreen et al.l 
I1992L Kendall et al. in preparation). The stellar mass 
dominates in the inner disk and gas falls into the po- 
tential well of the stellar spiral arms. These arms then 
propagate in the outer disk as sound waves. Since we 
reproduce the amplitudes of the stellar spiral wave cor- 
rectly, it is reasonable to believe that we are accurately 
representing the gas perturbations as well. 

To explore how the star formation in outer disks is 
affected by variations in spiral structure, we re-run our 
models again after increasing the halo mass by factors of 
two and four. As the halo to disk mass ratio increases, 
the disk becomes less self-gravitating and the power in 
the m = 2 mode is decreased, making the spiral struc- 
ture more flocculent. The structure of the disk is kept 
constant, and the number of particles in the halo is in- 
creased to keep the mass per particle in the halo constant 
between models. For the most massive disks, increasing 
the halo by four times means that the rotation curve 
peaks at over 300 km s _1 , which is unreasonably high, 
so these models are not included in the analysis. The 
m = 2 mode amplitude for the simulations with higher 
halo to disk mass ratios are shown on the middle and 
right panels in Figure [13] Note that as the halo to disk 
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Fig. 13. — The m=2 spiral mode amplitude for our models 
where the halo to disk mass ratio increases left to right. Disk mass 
fractions are 0.041 (left), 0.02 (middle) and 0.01 (right) of M to t. 
Colors and symbols correspond to those in Figure flTl 
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Fig. 14. — FUV images of the flat-mid run as the halo to disk 
mass ratio is increased by a factor of two (middle) and four (right). 
Contours are as in Figure [6] 

mass ratio increases, the m = 2 mode amplitude de- 
creases. T hese are still well within the constraints of ob- 
servations ([Rix fc Zaritskvlll995t lElmegreen et al.lll992l 
Kendall et al. in preparation). 

An example of the effect on the morphology of the 
outer disk star formation is shown in Figure [T4J Figure [14] 
shows the 2.2 Gyr snapshot for the flat-mid model as the 
halo to disk mass ratio is increased. The morphology 
is affected and this is reflected in the covering fraction. 
This indicates that in galaxies that are more dark matter 
dominated, such as dwarf galaxies, we might expect to 
see more flocculent spiral structures in the outer disk and 
smaller covering fractions. There is very little difference 
in other diagnostics. Plotting the star formation rate 
vs. gas surface density does not alter the results seen 
in Figure [TTJ In Figure [15] we show the T07 XUV type 
diagnostic diagram including the high halo mass models. 
The locus of points is almost entirely unaffected by the 
additional models. 

Changing the mass resolution of the halo particles also 
modifies the spiral structure. Higher resolution in the 
halo leads to less numerical noise and suppresses spi- 
ral structure. We choose our halo resolution because it 
recreates the observed amount of spiral structure. 

5.2. Stellar Profiles 

In a study of a lar ge sample of SDSS galaxies, 
iPohlen fc Truiillol (|2QQ6h found that only 10% of galax- 
ies have a purely exponential stellar profile. Most galax- 
ies show either a "down-bending break", where the ex- 
ponential profile of the outer disk is steeper than the 
inner disk, or an "up-bending break" where the expo- 
nential profile of the outer disk is shallower than that 
of the inner disk. While several other methods of pop- 
ulating the outer dis ks of galaxies have be en proposed 
(jRoskar et al.l l2008bl ]al: lYounger et all 120071 ). older stel- 
lar populations accompanying the low levels of star for- 
mation in outer disks have been suggested as an ex- 
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Fig. 15. — The properties of the T0 7 LS B zone for all models, 
including high halo mass models (see § 15. Note that the locus 
of points is similar to Figure [9] Colors and symbols are the same 
as Figure [9] 



planation for up-bending breaks (Thilk er et al.l l2QQ7t 
lElmegreen fc Hunter! |2QQ6[ ). The stellar profiles of our 
models, shown in Appendix A, exhibit features owing to 
over-densities, but generally appear to be characterized 
by a single exponential. The overall level of star forma- 
tion in outer disks is low, as shown by the star formation 
rate profiles in Appendix A, so it is not surprising that 
t he outer stellar profile chang es very little. 

lElme green fc Hunterl ([2006) present a model of star 
formation which creates down-bending breaks. These 
build stellar populations which also show a down-bending 
break. Their model includes a threshold for star forma- 
tion that depends on the Toomre stability criterion and 
a turbulent model of the ISM. Because our threshold for 
star formation is just a single volume density, the star 
formation profile of our simulations is determined en- 
tirely by the gas density profile of the galaxy. Our star 
formation rate profiles then flatten in the outer disk fol- 
lowing the flattening extended disk, rather than falling 
off with a sharper exponential. A more complex imple- 
mentation of the star formation threshold is needed in 
order to determine whether our models agree. 

5.3. Ionized Gas 

We construct the gas content of our model galaxies to 
reproduce observed HI distributions of galaxies. How- 
ever, there may be a significant amount of diffuse ionized 
gas, particularly in the outer regions of galaxies, that is 
not captured in 21 cm observations that measure the HI 
content. In the Milky Way the column densi ty ranges 
from 20% -60% of the column density of H I (|Revnoldsl 
Il99ll ). However, in outer disks, low gas densities reduce 
the efficiency with which the g as is able to shield itself 
from the UV bac kground (e.g. lElmegreen fc Parravanol 
119941 : ISchavell2004D . Studies of H I line- widths with radius 
indicate that t here is an unexplained source of heating 
in outer disks ([Tamburro et al.l [2009) which could also 
point to a contribution from the UV background. This 
indicates that H II could become a significant component 
of the ISM at large radii. If star formation truly does re- 
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Fig. 16. — Fraction of initial gas still present in outer disk as 
a function of time for the four simulations showing Type I XUV 
emission. Top Row: Mestel-mid and Mestel-high. Bottom: Flat- 
low and flat- mid. 

fleet the total gas mass, as it would if it depended on the 
mid-plane pressure of the ISM, the pressure owing to the 
HII mass could boost star formation rates in the outer 
disk. In the case of our simulations, adding HII would 
increase the amount of gas above the threshold in the 
outer disk and result in more star formation. Detailed 
studies of the ISM in outer disks are needed to determine 
whether this effect is important. 

5.4. Longevity 

There is much debate over whether star formation in 
outer disks is continuous or episodic. In many examples 
there is enough HI gas to support continuous star for- 
mation at the levels see n over the lifetime of the galaxy 
(|Gil de Paz et al.ll2QQ7h . Our simulations are consistent 
with this. In Figure [16] we show the fraction of the initial 
gas in the outer disk still present as a function of time for 
the four simulations that are identified as Type I XUV 
disks in § 14.11 The outer disk is defined to start where the 
star formation rate surface density falls below 3 x 10 -4 
M kpc -2 yr _1 in the initial disk. The oscillations owe to 
non-steady behavior in the disk as the simulation relaxes 
(see Figure [3]). Even the simulations with the most outer 
disk star formation, flat-mid and Mestel-high, lose only 
^ 20% of their gas over the 4 Gyr timeframe of the simu- 
lations. Flat-low and Mestel-mid show almost no change 
over this time. This indicates that XUV disk emission 
could continue for the lifetimes of the galaxies, though as 
the gas is consumed in the galaxies with higher levels of 
outer disk star formation, the amount of XUV emission 
may decline to what is seen in lower gas content models. 

5.5. Interactions 

In this paper, we study only isolated galaxies, where 
density perturbations in the outer disk result from spi- 
ral structure induced by numerical noise in the halo. Of 
course, other sources, such as interactions, can cause per- 
turbations in outer disks leading to in situ star formation. 
Strongly shocked gas owing to interactions in the outer 
disk would likely create high levels of star formation in 
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very small regions. We will explore this possibility in due 
course. 

5.6. Ha Emission 

While exploring the simulated Ha colors and distribu- 
tions of individual star clusters would be very interesting, 
our resolution limits us to particle masses of ~ 10 5 M . 
Since Ha knots in the outer disk are often illuminated by 
a single star emitting ionizing radiation, we do not have 
the resolution to study the Ha properties of our sources. 

6. CONCLUSIONS 

We have explored how outer disk star formation in 
extended gas disks resulting from a Kennicutt-Schmidt 
star formation prescription and volume density threshold 
compares to observations of XUV disks. We used a set of 
models with outer disk structures motivated by HI sur- 
veys. We find that models with slowly varying high sur- 
face density outer gas disks would be classified as Type 
I XUV disks. Possible variations in the star formation 
rate normalization or spiral structure amplitude mean 
that gas distributions that result in only a few knots of 
outer disk UV emission in our models may not always 
have UV emitting sources in their observed outer disks. 
While we cannot rule out other explanations for Type I 
XUV disks, since the Type I XUV emission in our models 
results naturally from empirical star formation laws and 
the observed structure of galaxy disks, we argue that the 
theory outlined here is a viable explanation. 

The star formation in the outer disk of our models is 
not prevalent enough to erase the tell-tale truncation of 
the profile caused by the star formation volume density 
threshold. Since B07 see no truncation in the UV with 
gas radial profile, this could indicate that we need to re- 
vise our star formation prescription. However, the star 
formation rates in our models are consistent with star for- 
matio n rates in the Type I XUV disk M 83 iDong et al.l 
(2008|. More detailed observational studies of the rela- 
tion of gas density to star formation rate in outer disks 
are needed to determine whether the observed fields of 
M 83 are a special case or whether disks selected for XUV 
emission are subject to different processes than the B07 
sample. 

While XUV emission does not depend enough on the 
detail of gas density profiles to distinguish easily between 
them, it is a good indicator of the presence of gas at 
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large radii in galaxies. Our results imply that nearly 
1/3 of galaxies have large amounts of gas at 2-3 times 
their optical extents. In addition, forming Type II XUV 
disks with traditional star formation laws must require 
some process, such as rapid gas accretion, for creating 
a large young stellar population without an accompany- 
ing older population. These results may be outlining a 
picture where accreting gas in the outer parts of disks is 
common at late times, as suggested by the inside out disk 
formation scenario d emonstrated by simulations of struc - 
ture formation (e.g. iMo et all 11993 iBrook et al1l2QQ6[ ). 
and supported by the possibility that cold clouds devel- 
oping in cosm ological simulations may accrete onto disks 
at late times ([Keres fc Hernquistll2QQ9[ ). 

The next step in exploring the observed signatures of 
in situ star formation would be to run very high resolu- 
tion simulations that could be compared to obse rvational 
work such as that of iHerbert-Fort et al.l (|2QQ9[ ). These 
authors have computed the cross correlation of young 
clusters in NGC 3184 and found a significant inter-cluster 
signal for separations of less than 1 kpc, with cluster 
masses of ~ 10 3 — 10 4 M . With higher resolution simu- 
lations, it is possible that measurements like these could 
differentiate between in situ star formation and stellar 
migration scenarios. Moreover, simulations with star 
formation prescriptions that are linked directly to the 
molecular gas may help r econcile in situ star formation 
and the results of B07. iRobertson fc Kravtsovl (|2QQ8[ ) 
present models indicating that directly relating the star 
formation rate to the molecular gas may eliminate trun- 
cations in the star formation rate profile. Continued 
observational studies and sophisticated simulations are 
needed to explore this complex regime in star formation. 
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Fig. 17. — Maps of (from top to bottom) stellar mass surface density, gas surface density, star formation rate surface density and plots 
of stellar profile (red) and gas profile (blue), star formation rate profile and Kennicutt-Schmidt law at 2.2 Gyr. The stellar mass profile 
from the initial simulation is over-plotted on the mass profile plot as the black dashed line. This is to highlight any differences from the 
original profile. The simulations are (left to right) exp-low with the disk mass equaling 0.041, 0.02 and 0.01 of the total mass, and exp-mid 
with a disk mass equaling 0.041 of the total mass. 
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Fig. 18. — Maps of (from top to bottom) stellar mass surface density, gas surface density, star formation rate surface density and plots 
of stellar profile (red) and gas profile (blue), star formation rate profile and Kennicutt-Schmidt law at 2.2 Gyr. The stellar mass profile 
from the initial simulation is over-plotted on the mass profile plot as the black dashed line. This is to highlight any differences from the 
original profile. The simulations are (left to right) exp-mid with the disk mass equaling 0.02 and 0.01 of the total mass, and exp-high with 
a disk mass equaling 0.041 and 0.02 of the total mass. 



17 




0.0 0.6 1.0 1.5 0.0 0-3 1.0 1.5 0.0 0.5 1.0 1.5 0.0 0.5 L0 1.5 
log 2 (W, pa"') log E (H» po~") log E flfo po"") log E (U B pa"") 



Fig. 19. — Maps of (from top to bottom) stellar mass surface density, gas surface density, star formation rate surface density and plots of 
stellar profile (red) and gas profile (blue), star formation rate profile and Kennicutt-Schmidt law at 2.2 Gyr. The stellar mass profile from 
the initial simulation is over-plotted on the mass profile plot as the black dashed line. This is to highlight any differences from the original 
profile. The simulations are (left to right) Mestel-low with the disk mass equaling 0.041, 0.02 and 0.01 of the total mass, and Mestel-mid 
with a disk mass equaling 0.041 of the total mass. 
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Fig. 20. — Maps of (from top to bottom) stellar mass surface density, gas surface density, star formation rate surface density and plots of 
stellar profile (red) and gas profile (blue), star formation rate profile and Kennicutt-Schmidt law at 2.2 Gyr. The stellar mass profile from 
the initial simulation is over-plotted on the mass profile plot as the black dashed line. This is to highlight any differences from the original 
profile. The simulations are (left to right) Mestel-mid with the disk mass equaling 0.02 and 0.01 of the total mass, and Mestel-high with a 
disk mass equaling 0.041 and 0.02 of the total mass. 
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Fig. 21. — Maps of (from top to bottom) stellar mass surface density, gas surface density, star formation rate surface density and plots 
of stellar profile (red) and gas profile (blue), star formation rate profile and Kennicutt-Schmidt law at 2.2 Gyr. The stellar mass profile 
from the initial simulation is over-plotted on the mass profile plot as the black dashed line. This is to highlight any differences from the 
original profile. The simulations are (Left to right) flat-low with the disk mass equaling 0.041, 0.02 and 0.01 of the total mass, and flat-mid 
with a disk mass equaling 0.041 of the total mass. 
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Fig. 22. — Maps of (from top to bottom) stellar mass surface density, gas surface density, star formation rate surface density and plots 
of stellar profile (red) and gas profile (blue), star formation rate profile and Kennicutt-Schmidt law at 2.2 Gyr. The stellar mass profile 
from the initial simulation is over-plotted on the mass profile plot as the black dashed line. This is to highlight any differences from the 
original profile. The simulations are (left to right) flat-mid disk with the disk mass equaling 0.02 and 0.01 of the total mass. 



